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Abstract 

Literature suggests that disease severity and neurotransmitter signaling patliway genes can accurately identify 
antipsychotic response in schizophrenia patients. However, putative role of signaling molecules has not been tested in 
schizophrenia patients based on severity of illness, despite its biological plausibility. In the present study we investigated 
the possible association of polymorphisms from five candidate genes RGS4, SLC6A3, PIP4K2A, BDNF, PI4KA with response to 
antipsychotic in variably ill schizophrenia patients. Thus in present study, a total 53 SNPs on the basis of previous reports 
and functional grounds were examined for their association with antipsychotic response in 423 schizophrenia patients 
segregated into low and high severity groups. Additionally, haplotype, diplotype, multivariate logistic regression and 
multifactor-dimensionality reduction (MDR) analyses were performed. Furthermore, observed associations were 
investigated in atypical monotherapy (n = 355) and risperidone {n = 260) treated subgroups. All associations were estimated 
as odds ratio (OR) and 95% confidence interval (CI) and test for multiple corrections was applied. Single locus analysis 
showed significant association of nine variants from SLC6A3, PIP4K2A and BDNF genes with incomplete antipsychotic 
response in schizophrenia patients with high severity. We identified significant association of six marker diplotype ATTGCT/ 
ATTGCT (rs746203-rsl0828317-rs7094131-rs2296624-rsll013052-rs1409396) of PIP4K2A gene in incomplete responders 
(corrected p-value = 0.001; adjusted-OR = 3.19, 95%-CI = 1 .46-6.98) with high severity. These associations were further 
observed in atypical monotherapy and risperidone sub-groups. MDR approach identified gene-gene interaction among 
eDA/F_rs7103411-SDNF_rsl491851-S/.C6/\3_rs40184 in severely ill incomplete responders (OR = 7.91, 95%-CI = 4.08-1 5.36). 
While /?GS4_rs2842026-S/.C6/l3_rs2975226 interacted synergistically in incomplete responders with low severity (OR = 4.09, 
95%-CI = 2.09-8.02). Our findings provide strong evidence that diplotype ATTGCT/ATTGCT of PIP4K2A gene conferred 
approximately three-times higher incomplete responsiveness towards antipsychotics in severely ill patients. These results 
are consistent with the known role of phosphatidyl-inositol-signaling elements in antipsychotic action and outcome. 
Findings have implication for future molecular genetic studies as well as personalized medicine. However more work is 
warranted to elucidate underlying causal biological pathway. 
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Introduction 

Schizophrenia is a heterogeneous mental illness with varying 
degrees of positive, negative and cognitive symptoms governed by 
both genetic and environmental factors. Accumulating literature 
suggests contribution of genetic variations may lead to differential 
antipsychotic response ranging from adequate response, to 
development of side-effects or treatment resistance [1]. Addition- 
ally affected individuals may have different combinations of 
disease symptoms which cause considerable disease heterogeneity 
and differential therapeutic outcome [2] . Polygenic involvement of 
genetic factors in the schizophrenia pathophysiology hinders the 
identification and characterization of causal genes that could 
develop as therapeutic targets. Inspite of the increasing number of 



pharmacogenomic studies attempts to explicidy uncover the 
mechanism behind the antipsychotic response remains inadequate 
which urges further investigation of candidate pathway genes 
including drug targets and signaling molecules since effect of 
antipsychotic drugs propagated through various biochemical 
processes that are part of coherent signaling network. Thus, it is 
important to investigate the influence of genetic variation and 
clinical factors on therapeutic outcome in differentially affected 
schizophrenia individuals according to the disease severity [3]. 
These facts prompted us to identify variants that may influence 
drug response by investigating critical signaling molecules known to 
express in different brain regions, implicated in schizophrenia and 
drug response. Among these, we prioritized Brain-derived neuro- 
trophic factor {BDMFj, Solute carrier family 6 (neurotransmitter 
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transporter), member 3 {SLCJiA'J), Regulator of G-protein signaling 
4 (RGS4), Phosphatidylinositol-5-phosphate 4-kinase, type II, alpha 
{PIP4K2A) and Phosphatidylinositol 4-kinase, catalytic, alpha 
[PMKA). 

The BDMF is known to influence neurogenesis, neuroplasticity 
and interact with other neurotransmitters like dopamine, gluta- 
mate, serotonin and gamma-aminobutyric acid. Polymorphisms in 
this gene have been studied for roles in efficacy and adverse effects 
of antipsychotics treatment [4]. Recently it is suggested that 
antidepressant and antipsychotic treatments ameliorates disturbed 
monoamine syst(;ms caused by reduced BDMF activity, while 
reduced BDNF levels lead to diflFerential treatment response in 
schizophrenia patients caused by stress and lack of neuronal 
activity [5]. SLC6A3 plays a critical role in controlling dopamine 
transmission (spatial and temporal domains) through the accumu- 
lation of dopamine in extracellular space which is major site of 
action of psychostimulant drugs [6-8]. Sjoholm et al., reported 
increased number of SLC6A3 binding sites in the schizophrenic 
patients who were being treated with dopamine D(2)-receptor 
blocking antipsychotics. Thus altered SLC6A3 functioning may 
lead to inter-individual variability towards antipsychotic drug 
treatment [9]. Whereas reduced RGS expression may result in 
prolonged signal transduction, thereby leading to altered neuro- 
transmission of dopamine, serotonin or glutamate. In addition, 
studies showed RGS4 is extremely responsive to external stimuli 
and hence can modulate GPCR neurotransmitter receptors by 
antipsychotics [10,11]. PIP4K2A and PMKA have been suggested 
as putative susceptibility genes on both positional and functional 
grounds as these are part of phosphoinositide signaling pathway 
which is implicated in schizophrenia etiology and may modulate 
antipsychotic response [12-16]. PMKA located at 22ql 1.2 deletion 
region repeatedly linked to schizophrenia other psychiatric 
diseases. While PIP5K2A maps to 10pl2 region that has been 
implicated in schizophrenia and bipolar disorder. These two 
kinases involved in the biosynthesis of phosphatidylinositol-4,5- 
bisphosphate (PIP2) by phosphorylation of inositol ring at D4 
position. PIP2, secondary messenger, is one of the key elements of 
crossroads in phosphoinositide signaling and plays a role in 
membrane transduction of neurotransmitter signals as well as in 
intracellular signaling pathways that are known to be implicated in 
schizophrenia [17]. 

Studies across populations could help to identify true genetic 
association but population structure and divergent linkage 
disequilibrium (LD) pattern have impact on associations which 
has lead to inconsistent results. Also, selection of markers under 
investigation could sufficiently contribute to non-replicability of 
associations. In addition to this, single marker expose only the 
localized effect if individual SNP being analyzed may not be the 
causal variant or in LD with the functional variant [18,19]. Thus 
haplotype (multimarker) analyses may overcome diflFerences 
among LD structure across populations and could be more 
informative in terms of identification of causal marker for disease 
development and therapeutic outcome [20-24]. 

With this background, in the present study, we investigated 53 
polymorphisms from five genes {RGS4, SLC6A3, PIP4K2A, BDNF 
and PMKA) in 423 south Indian schizophrenia cases segregated 
into low and high severity of illness to assess their influence 
on antipsychotic response. We observed association of six 
marker diplotype ATTGCT/ATTGCT (rs746203-rs 108283 17- 
rs7094131-rs2296624-rsll013052-rsl409396) of PIP4K2A gene 
with an estimated odds ratio of OR=3.19 (95%-CI= 1.46-6.98) 
in severely ill schizophrenia patients who had inadequate response 
for antipsychotics irrespective of treatment in the study population. 
This further suggests that endophenotypically distinct subgroups of 



schizophrenia patients should considered separately for personal- 
ized therapeutic approaches. 

Methods 

Participant recruitment and clinical assessment 

Our study is a pharmacogenetic based observational follow-up 
study conducted in a natural setting. The study sample comprised 
of 482 unrelated schizophrenia patients of south Indian origin, 
who met DSM-IV criteria for schizophrenia diagnosis (Table 1). 
These patients were recruited from the clinical services of 
NIMHANS during the period of March 2007 to March 2010. 
Best-estimate diagnosis was established on the basis of structured 
interviews by experienced psychiatrists using Schedules for 
Clinical Assessment in Neuropsychiatry (SCAN), OPCRIT 3.1 
[2.'5,26]. Patients were excluded if they had 1) past history or family 
history of any mental illness or neurological disorder, 2) history of 
substance dependence or abuse, 3) history of head injury or 4) 
pregnancy. Detailed clinical and treatment history including age of 
onset (AOO), duration of illness (DOI), medication, relapse, non- 
compliance and change of treatment were also documented. 
Further, patients were evaluated by clinical global impressions 
(CGI) scale at the time of enrolment and each component of the 
CGI was rated separately to assess the patient's status. This was 
followed by three months during which patients were treated with 
antipsychotic drugs. Medication was as chosen by the clinicians 
and administered drugs included typical antipsychotics such as 
chlorpromazine (200-400 mg/ day), fluphenazine (20-50 mg/ 
day), and flupenthixol (20-40 mg/ day), and atypical antipsychot- 
ics such as clozapine (25-400 mg/day), risperidone (2-10 mg/ 
day), olanzapine (5-20 mg/day), ziprasidone (40-160 mg/day), 
quetiapine (1200 mg/day), aripiprazole (10-45 mg/day), amisul- 
pride (100-600 mg/day), and levosulpride (200 mg/day) or 
combination of two or more atypical/typical antipsychotic drugs. 
After three months of treatment, patients were invited again and 
re-evaluated with CGI for antipsychotic response. Total 423 
patients completed three months follow-up period, among whom 
355 (83.92%) were on atypical monotherapy, and 260 (61.45%) 
were on risperidone (Table 2). It is expected that severely iU 
patients would have poor treatment outcome, especially if they do 
not respond to several months of treatment or if they are not 
responding, their symptoms may become more severe. Taking this 
fact into consideration, patients with CGI-Severity (CGI-S) score 
S3 (not at all ill to mildly ill) were categorized into low severity 
group (LSG) and those with CGI-S S4 (moderately to most 
extremely HI) were categorized into high severity group (HSG) [3] . 
Further, patients with CGI-Improvement (CGI-I) £2 or a drop of 
2 from baseline score were described as complete responders (CR) 
towards antipsychotic treatment. Patients with CGI-I S3 or more 
were grouped as incomplete responders (IR) [27]. Based on these 
criteria, 193 (45.63%) patients were grouped into LSG, among 
them 142 (33.57%) were CR and 51(12.06%) were IR and, 230 
(54.37%) constituted HSG group widi 68 (16.07%) CR and 162 
(38.30%) IR. Additionally 230 healthy individuals without any 
history of psychiatric illness were recruited from the same 
geographical area. Ethnicity of participants was assessed by self- 
report. 

Ethics Statement 

Written informed consent was obtained from all the participants 
after clinicians explained the study protocol to them. The study 
was approved by institutional review board of both the partic- 
ipating institutes: CSIR-Institute of genomics and integrative 
biology (CSIR-IGIB) and National Institute of Mental Health and 
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Table 1. Demographic and clinical characteristics of schizophrenia cases (n = 
study. 


482) and healthy controls (n 


= 230) enrolled in the 




Characteristics 
(mean±SD) 


Cases(n = 482) 


IVIales (n = 289; 60%) 


Female (n = 1 93; 40%) 


"P-value 


Age (yrs) 


29.52 ±7.33 


29.18±6.82 


30.06±8.05 


0.39 


Age at onset (yrs) 


24.98 ±7. 11 


24.73±6.41 


25.36±8.06 


0.96 


Duration of Illness (yrs) 


4.33±3.25 


4.33±3.52 


4.35 ±2.78 


0.51 


Disease Severity (CGI-S scale) 


3.53 ±1.43 


3.59±1.41 


3. 43 ±1.48 


0.23 




Controls (n = 230) 


Males (n = 142; 62%) 


Females (n = 86; 38%) 


"P-value 


Age (mean±SD) yrs 


29.48 ±8.51 


30.64±7.31 


27.61 ±9.01 


0.01 



Footnote: SD, Standard deviation. 

^P-values were calculated by using two-tailed Student's f-test. 
doi:l 0.1 371 /journal.pone.Ol 02556.t001 



Neuro Sciences (NIMHANS). Informed consent was obtained 
from respective family members if patient was not judged to be 
take liis or her own decisions or think appropriately. This 
assessment of clinical judgment was performed by primary care 
physicians and senior psychiatric clinicians. 

Population stratification 

A major concern of association studies was genetic heterogene- 
ity which may lead to false positive results [28] . Due to this reason, 
study participants were recruited from the same geographic 
region and were also genotyped for 10 autosomal microsatellite 
markers unlinked to schizophrenia for genetic homogeneity testing 
(Applied Biosystems linkage mapping set, version 2.0). These 
microsatellite markers were analyzed using GeneScan module 
of the Genotyper software, version 3.7 (Applied Biosystems). 
Stratification was tested by comparing genotype frequencies of 
each of the marker between patient and control groups using 
Pearson's tests. The observed test statistics p-value was 
estimated as the fraction of 10,000 simulated test statistics which 
exceeded the observed value. Further, for each locus, the sum of 
test statistics was computed with the number of degrees of freedom 
(df) which was equal to the sum of the number of df for the 
individual loci. This was performed by Structured Population 
Association Test (STRAT) program (Table Sla in File SI) [29]. 
In addition to unlinked microsatellite markers, we selected 441 
neutral markers to test presence of substructures in the studied 
population (Table Sib in File SI). These markers were selected 
from 4991 neutral markers from AfFymetrix 50K dataset [30,31]. 
Neutral markers which are not located in region linked to 
psychiatric illness and in genes known to be implicated in brain 
diseases, brain development and functioning were selected for 
present study. To test population stratification if any, the 
STRUCTURE software was used and admixture model was 
assumed [32]. According to this model individual i have inherited 
some fraction of his/her genome from ancestors in population K 
and it computes the proportion of the genome of an individual 
originating from each inferred population posterior mean 
estimates known as Q. This model is recommended due to its 
flexibility for dealing with real population complexities. We 
assumed K= 2 and model parameters includes length of burning 
period, number of Markov Chain Monte Carlo (MCMC) reps 
after burning, and iterations as 100, 100, and 100 respectively 
were considered. 



Selection of SNPs 

Most obvious candidates for pharmacogenetic studies are the 
genes that encode the drug targets. However, genes that might 
modulate drug-target interactions and downstream signaling 
cascade are also important [33]. Therefore, genes that are 
involved in neurotransmitter signaling and are known to modulate 
antipsychotic response namely RGS4, SLC6A3, PIP4K2A, BDNF 
and PI4KA were selected for present study. All the investigated 
polymorphisms were either previously reported and or were 
observed as functional in repositories such as, RegialomeDB and 
HaploReg that were used to screen the possible functional and 
regulatory effect of variants on gene regulation [34,35]. Non- 
coding SNPs located upstream and downstream of the gene 
spanning regulatory regions, predicted to be present in sphce site, 
and enhancer or silencer elements and untranslated region and 
coding SNPs (synonymous and non-synonymous) were selected 
to cover the entire gene. Total 1 1 7 SNPs were prioritized across 
five selected genes, of which 113 SNPs passed assay designing 
optimization and hence genotyped in entire sample set (Table S2 
in File SI). 

Genotyping 

Peripheral whole blood samples were obtained from all the 
participants by venepuncture, and genomic DNA was prepared 
using a modified salting-out procedure [36] . Prioritized SNPs were 
genotyped by primer extension reaction followed by the MALDI- 
TOF mass spectrometry (Sequenom). Automated calls were made 
using Mass Array Typer module of the software. To rule out any 
genotyping error, 5% of the total samples were re-genotyped 
randomly using single-base primer extension methods (SNaPshot, 
Applied Biosystems). Calls were made after analyzing data in 
PeakScanner Software (Applied Biosystems). Allele calls, clustering 
patterns and peak intensities were reviewed by an experienced 
technician. The SNPs with call rate exceeding 98% were 
considered. 

Statistical analyses 

Our analyses were focused on comparing the response of 
antipsychotics in presence of polymorphisms and adjusted for 
variables including gender, age, AOO, DOI, treatment type. 

Genotype-Phenotype association analyses. Genotyping 
was completed on a total of 1 1 3 SNPs; 1 2 SNPs had missing 
genotype data for all participants. The remaining 101, were tested 
for HWE using x""test at 0.001 level of significance and allelic 
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frequencies of all the variants were calculated for control and 
patient sample pools using statistical modules implemented in 
PLINK software [37]. Eight SNPs failed the HWE test (P<0.001) 
in controls and 17 SNPs failed in the minor allele frequency 
(MAF<0.1) test. Linkage disequilibrium (LD) among SNPs was 
estimated by correlation coefficient (r^) values for genotype data of 
all markers across the selected genes in controls using Tagger 
algorithm in Haploview program, version 4.1. [38]. Further, 23 
SNPs were pruned out based on LD (r^>0.9) with other genotyped 
SNPs in the studied population and 53 SNPs were included in 
analyses (Table S2 and S3 in File SI). 

All genotype-phenotype association analyses were performed on 
53 SNPs. In order to reduce number of statistical tests, dominant 
and recessive model of inheritance were selected to compare their 
distribution in LSG and HSG for treatment response in 423 
patients. These genetic models of inheritance were selected on the 
basis of distribution of genotypes and alleles of each SNP. To 
assess the study-wise significance, 10,000 permutation test were 
performed to obtained empirical p-values for all SNPs with the 
aim of calculating a value that controls for all the applied tests. 
Power of the study was calculated by PS Power and Sample Size 
Calculations, version 3.0, considering most significant p-value in 
single marker genotype-phenotype association analysis [39] . Entire 
genotype data was phased before constructing the haplotypes 
using the software PHASE 2.1 based on Bayesian algorithm. 
Parameter values of 100 iterations, a thinning interval of 10, and a 
burn-in value of 100 in the MCMC simulations were used [40]. 
SNPs found to be associated with drug response were included 
for haplotype analysis, using PLINK software and highly 
conservative Bonferroni method of multiple test correction was 
applied (p = 0.0027 = 0.05/18 p<0.05/n, where n = total number 
of statistical tests). Haplotypes that remained significant after 
correction were further examined for diplotype (homozygous 
condition). Additionally, p-values of single marker and diplotype 
were adjusted for gender, age, AOO, and DOI using multivariate 
logistic regression (STATA 1 1 .0). 

Subgroup analysis. The SNP, haplotype and diplotype that 
were observed to be associated with drug response in entire 
schizophrenia samples were also examined separately in subgroups 
with 355 patients on atypical antipsychotic monotherapy and out 
of these 260 patients were on risperidone. 

Epistasis analyses. Considering the heterogeneous antipsy- 
chotic response in schizophrenia patients one can expect 
interaction among the SNPs. Therefore, all the 53 SNPs were 
evaluated for potential interactions using multifactor dimension- 
ality reduction (MDR) method [41]. In brief, MDR is a non- 
parametric method that defines high and low-risk groups by 
combining lower and higher predisposing genotypes, which are 
then evaluated for its ability to categorize and predict drug 
response status through 10-fold cross validation. The cross 
validation consistency (CVC) was calculated for each possible 
model. Dendrogram were applied to portray SNPs interactions of 
the best model. 

Results 

Demographic and clinical characteristics of participants 

Out of 482 recruited schizophrenia patients 423 subjects 
(87.8%) had completed three months of follow up (male:253; 
59.8%, female: 170; 40.2%). The majority of patients were 
diagnosed with paranoid schizophrenia (Table 1). The demo- 
graphic characteristics of 423 schizophrenia patients eligible for 
genotype-phenotype correlation are summarized in Table 2. 
There was no significant gender diflFerence among CR and IR 
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Table 6. Distribution of diplotype, ATTGCT/ATTGCT {PIP4K2A) and its association analysis with drug response in schizophrenia 
patients {n = 423), atypical nnonotherapy {n = 355) and risperidone group (n = 260). 





Groups SNPs 


Number (frequency) 


OR (95%CI) 


P-value 




CR 


IR 






Schizophrenia patients 


Low severity rs746203|rsl0828317|rs7094131 |rs2296624|rs1 1013052|rsl409396 


33 (0.67) 


16 (0.32) 


1.51 (0.69-3.22) 


0.252 


High severity rs746203|rsl 082831 7|rs70941 31 |rs2296624|rsl 101 3052|rsl409396 


9 (0.14) 


55 (0.85) 


3.37 (1.51-8.28) 


0.001" 


Atypical monotherapy group 


Low severity rs746203|rsl 082831 7|rs70941 31 |rs2296624|rs1 101 3052|rsl409396 


28 (0.68) 


13 (0.31) 


1.34 (0.57-3.04) 


0.458 


High severity rs746203|rsl0828317|rs7094131 |rs2296624|rs1 1013052|rsl409396 


7 (0.14) 


43 (0.86) 


3.6 (1.45-10.15) 


0.002'' 


Risperidone group 


Low severity rs746203|rsl 082831 7|rs70941 31 |rs2296624|rs1 101 3052|rsl409396 


16 (0.59) 


11 (0.40) 


1.83 (0.66-4.90) 


0.18 


High severity rs746203|rsl 082831 7|rs70941 31 |rs2296624|rs1 101 3052|rsl409396 


4 (0.10) 


35 (0.89) 


5.09 (1.62-21.01) 


0.002'' 



Footnote: Haplotype withstand Bonferonni correction was included for diplotype analysis. Diplotype with minimum frequency less than 0.10 were ignored. 
"P-values less than threshold value of a which is equal to 0.05/6(n) = 0.008, where n is number of test for diplotype association analysis. 
CR, Complete responders; IR, Incomplete responders; OR, Odds ratio; CL Confidence interval. 
doi:1 0.1 371 /journal.pone.Ol 02556.t006 



in both the severity groups. In addition, no significant differences 
were observed in the age, AOO and treatment type among all the 
phenotypic groups studied except for DOI among CR and IR 
patients with high severity (unadjusted p-value —0.02). Addition- 
ally, treatment type was same in both severity groups. 

Population stratification results of ten unlinked microsatellite 
markers suggested homogeneity among analyzed sample sets with 
p-value of 0.142 (x^= 109.79; df=95) (Table Sla in File SI). 
Further, STRUCTURE analysis of 441 neutral markers revealed 
both cases and controls have same probability to be member of 
K= 2 clusters. Each individual in the data set represented by 
single vertical line divided into 2 colored segments {K= 2) that 
represents estimated membership fraction of each individual in 
both the K inferred clusters. Proportion of colored segments 
shows estimates of which is posterior mean estimate of each 
individual originating from each inferred population (Figure SI). 
In the summary plot depicts no clear distinction between the two 
assumed clusters hence we can conclude that studied samples were 
taken from same population pool. These results were further 
validated using STRAT software which computes summation of 
X^ = 353.603 (434 df), Global P-value after 1000 simulations for 
each locus found to be 0.9981 further supports that our studied 
sample pool is homogenous (Table Sib in File SI). 

Drug response analyses 

In 423 patients, 9 SNPs showed statistical significant difference 
among CR and IR in LSG or HSG with p-value <0.05 based on 
die best model fit (dominant or recessive model) (Table 3). 
Furthermore, p-values were adjusted for age, gender, DOI, 
AOO and drugs. Polymorphism rs 1409396 {PIP4K2A) showed 
most significant association (corrected p-value — 0.02, adjusted 
OR = 0.32, 95%-CI = 0.16-0.64) among complete and incomplete 
responders in HSG and remained significant after correcting for 
multiple comparisons using 10,000 permutations. Hence, patients 
in HSG having TT genotype are 68% less likely to respond 
towards antipsychotic drugs. The power of association was 0.963 
with studied sample size, true OR for drug response in severely ill 
IR relative to CR was 0.32 and the associated with this test of this 
nuU hypothesis was 0.05. Same associations were detected in 355 
atypical monotherapy and 260 risperidone sub-group. As expect- 
ed, the influence of particular drug was not observed because 



PIP4K2A is element of neurotransmitter signaling pathway and 
involved in both schizophrenia manifestation and downstream 
signaling cascade activated after drug-receptor interactions. 

Out of nine polymorphisms, eight SNPs (two of SLC6A3 and six 
of PIP4K2A) with p<0.05 were selected for haplotype analyses 
to examine their role in antipsychotics response. Remaining 
polymorphism rs56164415 was from BDMF gene and hence 
haplotype could not be constructed. Information for identified 
haplotypes is provided in Table 4. Five haplotypes were 
observed in 423 schizophrenia patients, among them, six 
marker 'ATTGCT' haplotype (rs746203-rsl0828317-rs7094131- 
rs2296624-rsll013052-rsl409396) of PIP4K2A gene showed 
significant association between CR (0.44) and IR (0.64) in HSG 
(p-value<0.002, OR = 2.05, 95%-CI= 1.26-3.33). In addition, 
this haplotype was associated with drug response in atypical 
monotherapy and risperidone subgroups and the association 
remained significant after applying Bonferroni corrections 
(Table 5). Further, this haplotype was examined in diplotype 
combination (ATTGCT/ ATTGCT) and found to be significantiy 
associated with incomplete antipsychotic response in patients 
with high severity (corrected p-value <0.0027; OR 3.37, 95%- 
CI= 1.51-8.28) (Table 6). Additionally, the effect of this diplotype 
on antipsychotics response in LSG and HSG was adjusted for age, 
sex, DOI and AOO in all patients, as well as in the treatment 
stratified groups using multivariate logistic regression analyses 
to negate any confounding factors. The analyses indicated 
independent effect of diplotype with drug response in severely 
iU patients among 423 samples (adjusted p-value = 0.004, 
OR=3.19, 95%-CI= 1.46-6.98), 355 atypical monotherapy 
(adjusted p-value = 0.004, OR = 3.66, 95%-CI= 1.51-8.88) and 
260 risperidone subgroup (adjusted p-value = 0.007, OR = 4.75, 
95 % -CI = 1.53-14.72). Receiver operating characteristic curve 
(ROC) clearly demonstrated that area under the curve was 
6% more in HSG (0.6, 95%-CI = 0.54-0.65) tiian LSG (0.54; 
95%-CI = 0.46-0.61) (Table 7). These results suggested an impact 
of this diplotype was larger in high severity patient than in patients 
with low severity indicating its role in therapeutic outcome. 

Epistasis analyses. The MDR analysis was performed in 
aU 423 patients to further investigate gene-gene interactions 
between CR and IR for both severity groups. Overall best models 
were iJr;,W_rs2842026-6ZC'&45_rs2975226 and 5i)7VF_rs7 10341 1- 
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Figure 1. Interaction dendrogram generated from analysis of 
53 polymorphisms by MDR method in patients with (A) low 
severity and (B) high severity of illness. A red or orange line 
indicates synergistic relationship; golden line represents additivity and 
blue or green shows redundancy. Short length arm of dendrogram 
indicates strong interaction. Footnote: IVIDR-Multifactor dimensionality 
reduction. 

doi:1 0.1 371/journal.pone.01 02556.g001 

BDJVF_rsl49mi-SLC6A3_rsWm with CVC of 7/10 and p- 
value<0.0001 in LSG and HSG patients respectively (Table 8). It 
was observed that in LSG, patients with AA genotype of rs2842026 
in combination with AA and AT genotypes of rs2975226 conferred 
4.09 times higher risk for responsiveness towards antipsychotic 
treatment (OR = 4.09, 95%-CI = 2.09-8.02), whUe HSG patients, 
AG and GG genotype of rs40 1 84 in combination with CT and TT 
genotype of rs7103411 and CT and TT genotypes of rsl491851 
posed lower risk for non responsiveness towards antipsychotic 
response (OR=7.91, 95%-CI = 4.08-15.36) (Fig. 1 and 2). 
These results support strong synergistic effect between these 
polymorphisms making an individual susceptible to incomplete 
antipsychotics response. 

Discussion 

Single nucleotide variations are known to account for 90% of 
genetic variability, hence, we examined their effect on antipsy- 
chotic response in schizophrenia patients with variable degree 
of severity. Among 53 SNPs investigated from five genes, we 
observed TT genotype of rs 1409396 [PIP4K2A) was significandy 
associated with incomplete antipsychotic response in severely iU 
patients in 423 schizophrenia patients as well as in subgroups. This 
polymorphism is present in intron 4 at conserved domain for 
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Figure 2. Distribution of high risk (dark shaded) and low risk 
(light shaded) genotypes of best models among the studied 
polymorphisms. Left bars illustrates distribution of complete 
responders and right bars shows incomplete responders for each 
genotype combination in (A) low severity and (B) high severity of illness 
patients with schizophrenia. 
doi:1 0.1 371/journal.pone.01 02556.g002 
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Figure 3. Identification of regulatory potential of six marker haplotype using ENCODE data from UCSC genome browser. 
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binding of CLOCK:BMAL transcription factor and 'T' allele is 
predicted to alter its binding site. This transcription factor widely 
participates in responses to environmental contaminants, hypoxia, 
neurogenesis, synaptic plasticity and circadian regulation [42]. In 
addition, nominal association was detected in major alleles of five 
variants from PIP4K2A (rsl0828317, rsl 1013052, rs2296624, 
rs7094131, rs746203) in severely iU patients who had incomplete 
response towards antipsychotic drugs. Furthermore, diplotype 
with major allele ^A-T-T-G-C-T' of six variants from PIP4K2A 
gene (rs746203-rs 1 08283 1 7-rs7094 1 3 1 -rs2296624-rs 11013052- 
rs 1409396) was found to be significantly associated which 
suggested that it may influence treatment non-responsiveness 
however, its actual elfects should be elucidated by in-vitro and 
in-vivo approaches. As predicted by ENCODE data (Regialo- 
meDb and HaploReg) variants constituting this diplotype have 
putative transcription factor binding sites which emphasizes that 
diplotype may have important biological consequences (Figure 3). 
Polymorphism rsl 1013052 was predicted to alter the motifs of 
ERalpha-a_disc4, Pbx3_disc2 which may influence the expression 
of estrogen alpha receptor and pre-B-ceU leukemia homeobox 3. 
The rs7094131 was found to be present in conserved region of 
intron 6 and may alter binding sites for various transcription 
factors whereas rs746203 was present in expression quantitative 
trait loci and DNAsel hypersensitive site. One non-synonymous 
variation rsl 0828 17 was known to down-regulate the synthesis of 
PIP2 in Xenopus system and completely abolish the kinase activity 
[43]. However, recently, investigators did not find significant 
difference in mutated-enzyme activity with wild-type enzyme [16]. 
Over-representation of this diplotype in severely ill incomplete 
responders is consistent with consequence of altered phosphoino- 
sitide pathway in schizophrenia symptomatology and drug 
response through modulation of genes comes to play after drug- 
target interaction. 

The PIP4K2A catalyses phosphorylation of PI(5)P (phosphati- 
dyUnosito-phosphate) and lead to formation of PI(4,5)P2 (PIP2) 
which further hydrolyzed by the action of phospholipase C 
into second messenger inositoltrisphosphate and diacylglycerol 
[44-46]. PIP2 was shown to be involved in the regulation of 
secretion, cell proliferation, differentiation, motility and, synaptic 
vesicle cycle [47]. It is a key player of membrane transduction of 
neurotransmitter signals, intracellular signaling, which may be 
impaired in schizophrenia. This pathway was proven to be 



potential target for lithium action for treatment of bipolar disorder 
[48]. PIP4K2A controls the functioning of KCNQ. channels that 
further suppress the basal activity of dopaminergic neurons 
and their firing in substantia nigra and ventral tegmental area of 
brain [43]. Moreover, PIP4K2A stimulate excitatory amino acid 
transporter, EAAT3 which clears glutamate from synaptic cleft. 
Impaired activity of PIP4K2A decreases EAAT3 activity and is 
expected to foster glutamate toxicity in schizophrenia patients 
[49]. 

To the best of our knowledge, this is the first study that 
examined PIP4K2A for antipsychotic response. Interestingly, we 
found similar frequencies of PIP4K2A variants upon comparison 
among all schizophrenia cases, patients in LSG and HSG with 
controls. We failed to detect association oi PIP4K2A in case-control 
analyses, irrespective of disease severity (data not shown). Our 
study further corroborated to existing evidence for involvement of 
PIP4K2A gene in schizophrenia disease and it was in conformity 
with other association studies where, major alleles of variants were 
found to be over-represented in disease individuals as compared 
to control [17,50,51]. Phased haplotype and diplotype are more 
informative and have high statistical power for an association than 
single SNP association because SNP do not arise independently; 
they have intrinsic dependency with one another due to their 
introduction in genealogical process [52]. Twelve locus haplotype 
oiPIP4K2A 'G-A-T-A-A-A-G-DEL-C-G-A-G' (rs7081744-rs943190- 
rs 1 08283 1 7-rs 1 08283 1 6-rs 1 0430590-rs746203-rs 1 053454-rs32 1 5826- 
rs8341-rs2364115-rs214836-rsl536332) showed over-transmission of 
major alleles to schizophrenia affected offspring of German 
and Israeli populations [17]. Similarly, major allele 'T' of 
rsl 08283 17 found to be associated with deficit schizophrenia 
and non-deficit schizophrenia patients in Dutch population [50]. 
Major allele 'C' of another intronic variant rsl 1013052 
was found to be over-transmitted in offspring affected with 
schizophrenia in Indonesian population [51]. On the contrary, a 
three-locus haplotype C-G-A (rsl 08283 17-rs746203-rs8341) of 
minor alleles was found to be associated in Chinese population 
[53]. Further, under-representation of common alleles preferen- 
tially in schizophrenia affected females was detected in Irish 
population [15]. The CC genotype (C is minor allele) of 
rsl 08283 17 was found to be associated in Russian population, 
however another study failed to detect association for PIP4K2A 
variant in German population [54,55]. These studies provide 
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additional evidence that defective phosphoinositide signaling 
contributes to schizophrenia etiology. We analyzed frequency of 
six variants of PIP4K2A gene across the globe and found similar 
frequency distribution as in study population (Table S4 in File 
SI). However, in African and Japanese population major allele 
of two variants (rs7094131; rs746203) were found to be under- 
represented. Intra-locus heterogeneity (different variants within 
the same gene) could be one of the reason for this variable allele 
frequency and PIP4K2A best explains this phenomenon. This 
heterogeneity has been proposed as a major fartor for the 
inconsistencies among the large population association test. 
Interestingly, these SNPs were not in LD with other remaining 
associated SNPs in 1000 genome populations as well as in our 
population. This indicated, that combination of these variants, if 
present may influence antipsychotic response in severely iU 
patients. Thus, we expect that, these polymorphisms may have 
role in gene expression and regulation and hence important for 
pharmacogenomic and other epidemiology studies. The largest 
prospective observational SOHO study suggested that it is 
important to assess treatment impact separately in dififerentially 
affected patients according to severity and type of schizophrenia 
symptoms [56]. Further, Chen and colleagues examined the 
magnitude of antipsychotic response to predict non-response as 
'not minimally improved', 'not much improved' and 'not 
remitted' in schizophrenia patients classified as moderately to 
severely ill and less than moderately Ul [3]. Our data further 
add to existing evidences indicating disease severity can 
modulate antipsychotic response. Thus, it would be helpful if 
analyses in future association studies consider clinical factors 
including disease severity separately for identification of predic- 
tive genetic markers that accurately defines treatment response. 

In addition to PIP4K2A, we also observed nominal significant 
association for rs56164415 (C270T, BDjXF), rs2550956 and rs6347 
(SLC6A3) with antipsychotics response. Although similar distribu- 
tion within other populations was observed for polymorphisms 
from RGS4 and PI4KA gene p-value <0.05 did not allow a 
genotype-phenotype correlation. This indicates that association 
diflference may be attributable to genetic heterogeneity among 
different ethnic groups. Inconsistent associations were reported 
for polymorphism rs56 16441,5 with antipsychotic response as it 
was not associated in Finnish population whereas haplotype 
230bp-C-A (GTn repeat-rs56164415-rs6265) was associated 
with risperidone response in Han-Chinese schizophrenia 
patients [57,58]. An updated meta-analysis revealed 'T' allele of 
rs56164415 associated with schizophrenia in East Asian popula- 
tion but not in Caucasian [59]. Also, rs56164415 predicted to be 
present in CpG island (high GC content-63.6%) can influence 
gene expression and post-transcriptional regulation of gene 
expression [60]. Other than this, rs6265 and rsl 1030104 were 
extensively studied in context of disease and antipsychotic r(^sponse 
but, so far associations are inconclusive [61,62]. Reducxxl BDJ\T 
levels were observed in medicated and drug naive schizophrenia 
patients as indicated by meta-analysis with unexplained population 
heterogeneity [63]. In a recent study, authors detected overall 
facilitation of Bdnf expression in aripiprazole-treated animals and 
this provide novel information regarding the mechanisms that 
aripiprazole may regulate brain function through BD}fF and could 
contribute to improvement of schizophrenia associated neuroplas- 
tic [64]. SLC6A3 variants rs2550956 and rs6347 were associated 
with schizophrenia [65-68] . These multiple genetic variants might 
have strong synergistic effects beyond their weak individual effect. 
This indicates that inappropriate dopamine efflux may be an 
unexplored mechanism that modulates the therapeutic actions of 
existing dopamine and serotonin receptor antagonist drugs. Even 



though some candidates from pharmacogenetic studies were 
replicated and seemed to be promising, still to date, there are no 
predictive marker developed that are commonly accepted to assess 
antipsychotic response. Heterogeneity among genetic architecture, 
populations and phenotype, different diagnostic criteria, gene-gene 
and gene-environmental interactions could be reasons for 
the failure of pharmacogenomic studies. In present study, 
MDR results suggested interaction between RGS4_rs2M2Q2&- 
SLC6A3_rs2975226 in schizophrenia patients of LSG and between 
BDNF_r%l 1 034 1 1 -BDNF_n 1 49 1 85 1 -SLCGAS^rsW 1 84 in severely 
ill patients for their inadequate treatment response. These 
epigenetic effects could best explain descriptive phenotyjies such 
as complete responders and incomplete responders, variable 
degree of severity. This also reflects pleiotropic effects of variants 
which may have important implications for drug discovery. 

Observed associations in the present study are true for several 
reasons as we applied stringent quality control procedures for 
inclusion of SNP, performed diplotype analysis with appropriate 
multiple test corrections and study participants were genotyped 
with microsateUite markers for testing homogeneity. However, 
results of the present study need to be interpreted cautiously 
because of three main limitations. First, limited number of 
candidate genes was analyzed. Secondly, our sample pool was 
in-sufficient for testing association of rare and uncommon variants. 
However, common variants have comparatively high power and 
accounts for variabiUty in endophenotype for a disease [69]. The 
third hmitation was lack of replication of results in another cohort. 
Studies with larger coverage and gene mapping with comprehen- 
sive analysis and replication are required to support and provide 
biological insight into disease pathogenesis and drug response 
mechanism. 

In conclusion, study highlights the importance of PIP4K2A as 

crucial signaling partner in drug induced signaling and its 
relationship with clinical characteristics of disease as we identified 
a risk diplotype of PIP4K2A which may predict antipsychotic 
incomplete response in severely iU schizophrenia patients. The 
identified marker may be used for early detection of patients who 
are not responding to their current treatment and may require 
change for optimum therapeutic regimen and can be exploited 
for novel therapeutics as add-on to antipsychotics for betterment of 
the patients. However, study merits replication of results in a 
larger population followed by functional validation. The complex 
etiology of schizophrenia merits the consideration of both genetic 
and epigenetic systems and the meticulous experiment designs that 
unravel the underlying mechanisms conferring accountability for 
schizophrenia and development of new efficacious therapies. 

Supporting Information 

Figure SI Summary plot of estimates of Q. Single vertical line 
represent each individual. line broken into K colored segments 
with lengths proportional to each of the K inferred clusters 
including data of 441 neutral markers. 

gpG) 

File SI Table Sla, Genetic homogeneity test for stratification 
by comparing distribution of 10 unlinked microsatellites markers 
between schizophrenia patients (n = 482) and healthy controls 
(n = 230). Footnote: x^-Chi Sqaure; df-Degree of freedom. Table 
Sib, Genetic homogeneity test for stratification by comparing 
distribution of 441 bi-alleUc neutrsd markers between schizophre- 
nia patients (n = 482) and healthy controls (n = 215). Table S2, 
Number statistics of prioritized SNPs for five genes. Footnote: 
SNPs- Single nucleotide polymorphism; HVVE- Hardy-Weinberg 
equilibrium; LD- Linkage disequilibrium; MAF- Minor allele 
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frequency; r^- correlation coefficient, measurement of LD. Table 
S3, Details and allele frequencies of polymorphisms in schizo- 
phrenia patients (n = 482) and controls (n - 230). Footnote: ^From 
RegulomeDB and HaploReg; ^'Major allcle/minor allele; 'Minor 
allele frequency. Table 84, Global frequency and world-wide 
frequency of major allele six drug response associated polymor- 
phisms of PIP4K2A typed in 1000 genomes. Footnote: N-number 
of individuals; Major allele and major allele frequencies were 
provided in this table. 
(DOC) 
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